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We present our analysis of the decay constants of the beauty vector mesons B* and B* within the 
framework of dispersive sum rules for the two-point correlator of vector currents in QCD. While the 
decay constants of the vector mesons /s* and fs* —similar to the decay constants of the pseudoscalar 
mesons fs and /b„— individually have large uncertainties induced by theory parameters not known 
with a satisfactory precision, these uncertainties almost entirely cancel out in the ratios of vector over 
pseudoscalar decay constants. These ratios may be thus predicted with very high accuracy due to the 
good control over the systematic uncertainties of the decay constants gained upon application of our 
hadron-parameter extraction algorithm. Our final results read /s* //s = 0.944±0.011ope ±0.018syst 
and /bj /fse = 0.947 ±0.023opb ±0.020syst. Thus, both /a* //s and fs* / fB„ are less than unity at 
2.5(t and 2a level, respectively. 

PACS numbers: 11.55.Hx, 12.38.Lg, 03.65.Ge 


1. INTRODUCTION 

The QCD sum-rule approach [H-Q, based on the application of Wilson’s operator product expansion (OPE) to the 
properties of individual hadrons, has been extensively used for predicting heavy-meson decay constants. An important 
finding of these analyses was the strong sensitivity of the decay constants to the values of the input OPE parameters 
and to the prescription of fixing the effective continuum threshold [^. The latter governs the accuracy of the quark- 
hadron duality approximation and, to a large extent, determines the extracted value of the decay constant. Even if the 
parameters of the truncated OPE are known with arbitrarily high precision, the decay constants may be predicted with 
only limited accuracy, which we refer to as their systematic uncertainty. In a series of papers @ , we have formulated 
a new algorithm for fixing the effective threshold within Borel QCD sum rules and for obtaining reliable estimates for 
the systematic uncertainties. This procedure opened the possibility to provide predictions for the decay constants with 
a controlled accuracy i0- 

Here, we study the decay constants of the vector beauty mesons fs* and fs* ■ As is already known from the analysis 
of the decay constants of the pseudoscalar mesons B and Bg Q , the OPE uncertainties in the obtained predictions are 
rather large. The same occurs also for the B* and B* mesons. However, the OPE uncertainties to a great extent cancel 
out in the ratios of the decay constants of vector and pseudoscalar beauty mesons. An important result reported here 
is that the systematic uncertainties of the decay constants are rather small and well under control. Therefore, these 
ratios are predicted with a very good accuracy. It should be taken into account that we address a rather subtle effect 
at a few-percent level; a priori, it is not clear whether QCD sum rules are, in principle, capable to provide theoretical 
predictions at this level of accuracy. Obviously, the control over the systematics is becoming crucial. 

The ratio of the decay constants of vector over pseudoscalar heavy mesons is an interesting quantity: it is known to 
be unity in the heavy-quark limit and to approach this limit from below because of the radiative corrections @ . For 
beauty mesons, the few existing sum-rule analyses (which, however, could not gain good control over the systematic 
uncertainties) reported /b* /Jb slightly above unity din. Constituent-quark models typically also yield Jb* j! b> ^ 
[HI- A similar conclusion has been reached by interpolation of the lattice data from the charm-quark mass region to 
the beauty-quark mass [l^ • 

The first indication that this ratio for beauty mesons is below unity was given in our paper s 11311 . Recently, HPQCD 
0 also reported an accurate value of Jb* I Jb < in excellent agreement with the results of [l3j| . The analysis of |13l |. 
although conclusively indicating Jb* //b < 1, observed an unpleasant dependence of the extracted decay constants of 
the vector beauty mesons on the renormalization scale fj, chosen for the evaluation of the vector correlation function. 
This analysis solves the problem of the sensitivity to the choice of the scale /i by improving the extraction procedures 
for the decay constants and arrives at new predictions stable with respect to the choice of fi. Our detailed results read 

/b*//b = 0.944 ±0.011opE±0.018,y,t, /b|//b. = 0.947 ± 0.023ope ± 0.020,y,t, (El) 

in more than excellent agreement with the latest results from lattice QCD 0 . Let us emphasize once more that the 
OPE uncertainties cancel to a large extent in the above ratios. Thus, decisive for obtaining an accurate sum-rule result 
is our capability to control the systematic uncertainties of the QCD sum-rule method. 
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2. QCD VECTOR CORRELATOR AND SUM RULE FOR VECTOR-MESON DECAY CONSTANT fv 


The decay constants of ground-state vector mesons may be extracted by analyzing the two-point correlation function 

* J d^xe^P^{0\T{j,{x)jU0)) | 0 ) = ( 2 . 1 ) 

of the heavy-light vector currents for a heavy quark Q of mass mq and a light quark q of mass m, 

j/iix) = q{x)j^Q{x), (2.2) 

or, more precisely, the Borel transform of its transverse structure n(p^) to the Borel variable r, n(T). Equating n(T) as 
calculated within QCD and the expression obtained by insertion of a complete set of hadron states yields the sum rule 

CO CO 

U{t) = fyMye~‘^^'" + J dse“®’‘phadr(s) = J ds e“®^Ppert(s, m) -b npower(r,/r). (2.3) 

Sphys {mQ+my 


Here, My labels the mass, fv the decay constant, and e^(j>) the polarization vector of the vector meson V under study: 

{0\ql^l.Q\V (p)) = fvMye^[p). (2.4) 

For the correlator (EH), Sphys = + is the physical continuum threshold, wherein Mp denotes the mass of the 

lightest pseudoscalar meson containing Q. For large values of r, the ground state dominates the correlator and thus 
its properties may be extracted from the correlation function EH- 

In perturbation theory, the correlation function is found as an expansion in powers of the strong coupling “constant” 
Qfs (m) . The best known three-loop perturbative spectral density has been calculated in [l^ in terms of the pole mass of 
the heavy quark Q (that is, in the present case, Mb) and for a massless second quark [hereafter, we use the abbreviation 
a(h') = as(j^)/7r, where as{p) is the running coupling at renormalization scale v in the MS scheme]: 

Ppert(s) = (s, Mb) + a{v)p^^^ (s, Mb) + a'^{v)p^‘^'' (s, Mb, m) H-. (2.5) 

For both quarks having nonzero masses, the two-loop spectral density in terms of their pole masses was obtained in Q. 

The power corrections are also separately scale-independent; their explicit expressions can be found in [^, [l^. For 
instance, for pseudoscalar (P) and vector (U) currents the quark-condensate contributions may be written in the form 


npower(T) = -mb{y){qq{v))Mb 


exp(-M2r) 1 + -Cpa - -CFar(0,M,V) 


hL 


At) = -mb[v){qq{v)) 


1 


1 


exp(—M^r) 1 -b -Cpa 1 -b -C'_Fa r r(— 1, M^r) 


( 2 . 6 ) 

(2.7) 


where mb{v) is the 6-quark MS mass at renormalization scale v, mb{v){qq{v)) is a scale-independent combination, and 
r(n,z) is the incomplete gamma function El- 

However, even if the lowest-order contributions to the perturbative expansion and the vacuum condensates of lowest 
dimensions are known to good accuracy, a truncated OPE does not allow one to calculate the correlator for sufficiently 
large r, such that the continuum states give a negligible contribution to n(r) in the corresponding range of r. In order 
to get rid of the continuum contribution, the concept of duality is invoked: Perturbative-QCD spectral density Ppert(s) 
and hadron spectral density Phadr(s) resemble each other at large values of s; thus, for values of the integration lower 
limit s chosen sufficiently large, that is to say, (far) above the resonance region, one arrives at the duality relation 

OO OO 

J dse“®^Phadr(s) = J dse“®^Ppert(s). (2.8) 

s s 

Now, in order to express the hadron continuum contribution in terms of the perturbative contribution, the relation (12.811 
should be extended down to the hadronic or physical threshold Sphys. However, since the spectral densities Ppert(s) and 
Phadr(s) obviously differ in the region near Sphys, one can reasonably only expect to obtain a relationship of the form 

OO OO 

j dse“®^/9hadr(s) = j rfs e“®^Ppert (s) , (2.9) 

-Sphya Seff(T) 
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where the effective threshold Seff (t) is clearly different from the physical threshold Sphys, Seff (t) ^ Sphys, and, moreover, 
must be a function of the Borel parameter r Si- By virtue of (12.91) . we may hence rewrite the QCD sum rule (12.3p as 

Settl-r) 

= j ds e“'’'^Ppert(s, + npower(r, = ndual(T, Seff (t)). (2-10) 

(mQ+m)^ 

We refer to the right-hand side of this relation as the dual correlator^ and to the masses and decay constants extracted 
from this expression as the corresponding dual quantities. In addition to Ppert(s, d) and npower('7‘) m); the extraction of 
fv requires, as further input, a criterion that fixes the functional behaviour of the effective continuum threshold Seff (t). 

We shall demonstrate that QCD sum rules allow a very satisfactory extraction of the vector-meson decay constants, 
with an accuracy that is certainly competitive to that found within the framework of lattice QCD. 


3. OPE AND CHOICE OF RENORMALIZATION SCHEME AND SCALE FOR HEAVY-QUARK MASS 

The starting point of our discussion is the OPE for the correlator m- The three-loop perturbative spectral density 
/Opert(s; M) was calculated in [l^ in terms of the pole mass of the heavy quark. A nice feature of the pole-mass OPE is 
that each of the known perturbative contributions to the dual correlator is positive. Unfortunately, the pole-mass OPE 
does not provide a visible hierarchy of the perturbative contributions to the extracted predictions, which raises doubts 
whether the 0(a^)-truncated pole-mass OPE is indeed a good starting point for a reliable analysis of decay constants. 

A well-known remedy is to reorganize the perturbative expansion in terms of the &-quark running MS mass 
related (in the notations of M) to the corresponding pole mass Mb by 

Mb = rfibi^i)/ (^1 -I- -I- -I- O(a^). (3.1) 

The spectral densities in the MS scheme are found by expanding the pole-mass spectral densities in powers of a(/i) and 
omitting terms of order 0{a^) and higher; starting at order 0{a), they contain two parts: the “genuine” part from [Tsi l 
and the part induced by the lower perturbative orders when expanding the pole mass in terms of the running mass. 
By this, however, dne to the truncation of the perturbative series, one gets an explicit (unphysical) dependence of the 
dual correlator and of the extracted decay constant on the scale In principle, any scale should be equivalently good. 
In practice, however, the distinctness of the hierarchy of the perturbative contributions to the dual correlator depends 
on the precise choice of the scale. This opens a possibility of choosing the scale fi such that the hierarchy of the new 
perturbative expansion is improved. 

Figures [Hand [2] depict the dual decay constants of the B* and B mesons, respectively. For the &-quark MS mass, 
we use the value determined in by matching our QCD sum-rule results for fs to those of lattice QCD:^ 

mb(mb) = (4.247 ± 0.034) GeV. (3.2) 

The numerical values adopted for other relevant OPE parameters are BE El, in 

to( 2 GeV) = (3.42 ± 0.09) MeV, TOs( 2 GeV) = (93.8 ± 2.4) MeV, a^iMz) = 0.1184 ± 0.0020, 

= (0.024 ±0.012) GeV'‘, (g9)(2 GeV) =-[(267 ± 17) MeV]^ = 0.8 ± 0.3. (3.3) 

\ TT / {qq)(2 GeV) 

The purpose of Figs. |T] and H] is the illustration of the main features of the dual correlators (12.101) . therefore the QGD 
sum-rule estimates shown here are obtained for a r-independent effective threshold: Seff = const. The numerical 
value of the latter is, in each case, found by requiring maximal stability of the extracted decay constant in the 
Borel window. We emphasize that our results for the decay constants reported in the next Sections are obtained using 
the r-dependent effective thresholds. 

From Figs, d] and m we conclude that the 0(Q;g)-truncated pole-mass OPE exhibits no hierarchy of the perturbative 
expansion and better should not be used. Unfortunately, the hierarchy of the running-mass OPE is also not guaranteed 
automatically and depends strongly on the scale /i. 


^ As shown in [3, the PGD average mb{mb) = (4.18±0.030) GeV (see also for a recent overview of the 6-quark mass results) leads 
to a considerably larger value of fs , incompatible with the latest lattice-QCD results. However, the precise value of m 5 (m 5 ) has negligible 
impact on the ratio of the decay constants of vector and pseudoscalar mesons. 
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Fig. 1: QCD sum-rule estimates of the 5*-meson decay constant using the pole-mass OPE (a) and the running-mass OPE at the 
renormalization scales /r = 2.5 GeV (b), /r = 3 GeV (c), and = 5 GeV (d). The running-mass OPE for rnb(rnb) = 4.247 GeV is 
shown. The pole-mass OPE employs the corresponding two-loop pole mass Mb = 4.87 GeV. For each case, separately, a constant 
effective continuum threshold SeU is determined by requiring maximal stability of the predicted decay constant in a Borel window 
of the maximal width 0.05 < r (GeV~^) < 0.15. Bold lines (lilac)—total findings, solid lines (black)—0(1) contributions; dashed 
lines (red)—0(as) contributions; dotted lines (blue)— 0{ag) contributions; dot-dashed lines (green)—power contributions. 
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Fig. 2: Same as Figure[T]but for the B meson. 
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Let us define a scale jl by demanding = rnb{jl). From the 0{a^) relation between MS and pole mass, we find 
ft Ki 2.23 GeV. At this scale, the perturbative hierarchy of the MS expansion is worse than that of the pole-mass 
expansion because the 0(1) spectral densities coincide, whereas the 0{as) spectral density in the MS scheme receives 
a positive contribution compared to the pole-mass scheme. For lower scales fj, < fi, the hierarchy of the MS expansion 
gets worse with decreasing /i. For higher scales fi > fi, first the hierarchy of the MS-expansion improves with rising /r 
(Figs. [Hand[2]) . However, as the scale n becomes sufficiently larger than jl, the “induced” contributions, which mainly 
reflect the bad-behaved expansion of the pole mass in terms of the running mass, start to dominate over the “genuine” 
contributions. This is evident in Figs. [T] and HJ at /x = 5 GeV, the 0(1) contribution to the dual correlator rises steeply 
with T, whereas the Oja) contribution becomes negative in order to compensate the rising 0(1) contribution. Finally, 
for large values of fi we mainly observe a compensation between the “induced” contributions. We may expect in this 
case the accuracy of the expansion to deteriorate. 

Figures [Hand ID also reveal an essential difference between pseudoscalar and vector correlators: at the same scale /x, 
the good reproduction of the observed mass of the vector meson requires lower values of r compared to its pseudoscalar 
partner. This implies that the Borel window for the vector correlator should be chosen at lower values of r than the 
corresponding window for the pseudoscalar correlator. Moreover, for /x > 5-6 GeV the vector-meson mass cannot be 
reproduced in a reasonably broad r window and so the QGD sum rule cannot predict the vector-meson decay constant. 

For the present analysis, we thus choose as range of scales /x = 3-5 GeV: On the one hand, in this range we observe 
a reasonable hierarchy of the perturbative contributions to the correlator. On the other hand, we shall see that for this 
range of scales one can find sufficiently broad r windows where the decay constants may be reliably extracted by our 
algorithm. For the vector mesons, the upper bound of this window depends on /x. 


4. EXTRACTION OF THE BEAUTY-MESON DECAY CONSTANTS FROM OUR QCD SUM RULE 


In order to extract the decay constants, we first have to find a r window such that the OPE provides a sufficiently 
accurate description of the exact correlator (i.e., all higher-order radiative and power corrections are under control). 
Next, we must determine the r dependence of the effective threshold Seff(T). The appropriate algorithm was developed 
and verified within quantum-mechanical potential models Q and shown to work successfully for the decay constants of 
heavy pseudoscalar mesons @ . We introduce a dual invariant mass Mduai and a dual decay constant /dual by defining 

^dual('r) = logHduaKr, Seff(r)), fLA^) = ^Hdua^T, Seff(r)). (4.1) 

For a properly constructed HduaK^j Seff(T)), the dual mass coincides with the actual ground-state mass My. Therefore, 
any deviation of the dual mass from My is an indication of the contamination of the dual correlator by excited states. 

For any trial function for the effective threshold, we derive a variational solution by minimizing the difference between 
the dual mass dm and the actual (i.e., experimentally measured) mass in the Borel window. This variational solution 
provides the decay constant then via dUD- We consider a set of polynomial Ansatze for the effective threshold, viz., 

n 

(“- 2 ) 

j=o 


and fix the coefficients 


(the knowledge of which then allows us to compute the decay constant fy) by minimizing 


= (4-3) 

i=l 

over the Borel window. Still, different Ansatze for Sefi(T) yield different sum-rule predictions for the decay constants. 

Gareful studies of quantum-mechanical potential models indicate that it suffices to allow for polynomials up to third 
order: In this case, the band delimited by the predictions arising from linear, quadratic, and cubic Ansatze for Sefi(T) 
encompasses the true value of the decay constant. Even a good knowledge of the truncated OPE does not allow us to 
determine the decay constant precisely but it enables us to provide a range of values containing the true value of this 
decay constant. The width of this range may then be regarded as the systematic error related to the principally limited 
accuracy of QCD sum rules. Presently, we are not aware of any other possibility to acquire a more reliable estimate for 
the systematic error. Noteworthy, considering a r-independent threshold would not allow us to probe the accuracy of 
the obtained estimate for fy. 

On top of the systematic error comes the OPE-related error of the decay constant: the OPE parameters are known 
only with some errors, which induce a corresponding error of fy. We determine this OPE-related (statistical) error by 
averaging the results for the decay constant assuming for the OPE parameters Gaussian distributions with the central 
values and standard deviations quoted in (13.31) and a flat distribution over the scale /x in the range 3 < /x (GeV) < 5. 
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A. Decay constant of the f?* meson 

1. Choice of renormalization scale 

In principle, the decay constant should be independent of the scale fj, at which the correlation function is evaluated. 
In practice, however, due to the truncations of the perturbative expansion and the series of power corrections, and the 
neccessity to isolate the ground-state contribution from the hadron continuum states, a reliable extraction of the decay 
constant may be performed in only a limited range of the scale /r. For the vector beauty meson, the suitable range of 
is found to be ^ = 3-5 GeV: For /r < 3 GeV, the perturbative expansion for the vector correlator does not exhibit a 
satisfactory perturbative convergence and therefore gives no reason to believe that the unknown higher-order radiative 
corrections both in the perturbative part of the correlation function and in the radiative corrections to the condensates 
are negligible. At higher scales /r > 5 GeV, the B* mass cannot be reproduced with the required accuracy, signalling 
that there the contamination of the excited states cannot be cleaned out. 


2. Choice of Borel-parameter window 

We require that the B*-B mass splitting and the masses of B* and B mesons are reproduced, separately, with an 
accuracy not worse than 5 MeV for any r value within the selected ranges. As follows from the properties of the dual 
correlators, this requirement provides two constraints on the choice of the r window for B*: 

1. The T window for B* should be chosen at lower values of r compared to the B-meson case. 

2. The precise choice of the r-window for B* should correlate with the scale /r at which the correlator is evaluated. 

To satisfy the above criteria for B*, we set the lower boundary at Tmin (GeV“^) = O.OI and choose a /r-dependent upper 
boundary of the form Tmax (GeV“^) = 0.31 — 0.05 ^ (GeV), which choice enables us to extract Jb* with a systematic 
uncertainty not worse than 5 MeV and strongly diminishes the unphysical scale dependence of the decay constant /b» . 

Figure [3] shows the application of our procedure for fixing the effective threshold and extracting the resulting Jb*- 
The dependence of our QGD sum-rule result on the relevant OPE parameters, i.e., the &-quark mass mt = rni,{rnb), the 
quark condensate (qq) = {qq{2 GeV)) and the gluon condensate (aGG), proves to be well described by a linear relation: 

{qq), {aGG)) = (181.8 ± fy,,t) (l + (l “ MeV, (4.4) 

with 

,, mb-4.247 GeV ^ _ \{qq)\^/^ - 0.267 GeV ^ _ (aGG) - 0.024 GeV^ 

“ 0.034 GeV ’ ~ 0.017 GeV ’ “ 0.012 GeV'‘ ' ’ 

The above parameters 6 take values between —1 and -|-1 when the corresponding OPE quantity varies in its Icr interval. 
Varying all other OPE parameters in their Ict ranges leads to an effect on /g. of less than 1 MeV and is not shown here. 
Trusting in our experience from exactly solvable examples, we assume that the systematic uncertainty interval contains 
the true value of the decay constant and that inside this interval the true decay-constant value has a flat distribution. 

As evident from Eig.[3Ka), using a constant threshold leads to a contamination of the dual correlator by excited states 
(beyond the acceptable level), while this contamination is strongly reduced for n > 0: the values of the decay constant 
in Eig.[3Kb) resulting for n > 0 are nicely grouped together, whereas the n = 0 prediction emerges some 10 MeV below. 

A particularly convincing feature of the presented extraction procedure is the insensitivity of the extracted value of 
/b* (as well as that of Jb) to scale variations in the interval /r = 3-5 GeV (Fig.0]), achieved by demanding an accurate 
reproduction of the B* mass in the full t window, which requires a specific choice of the r window correlated with the 
scale fjL at which the correlator is evaluated. Such choice of the r window allows us to keep the systematic uncertainty, 
estimated by the half width of the band encompassing the results for the linear, quadratic, and cubic thresholds, at a 
level below 4 MeV in the full ^ range. Therefore, (14.411 describes well the result for any /i from the range n = 3-5 GeV. 

Assuming Gaussian distributions for all OPE parameters collected in (13.31) . we get the distribution of /b* depicted in 
Fig. [S] For the average and the standard deviation of the i?*-meson decay constant, we obtain 

/b. = (181.8 ± 13.10PE ± 4,y,t) MeV. (4.6) 

The OPE uncertainty is composed as follows: 11 MeV are due to the variation of mb and 6 MeV arise from the quark 
condensate. The uncertainties of all other OPE parameters contribute less than 1 MeV to the OPE uncertainty of /b*. 
The corresponding QGD sum-rule outcome for the B-meson decay constant /b from our earlier investigation [l8l| reads 








Mdual/Mv 
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Fig. 3: Dependence on the Borel parameter r of the dual mass (a) and the dual decay constant (b) of the B* meson, obtained by 
adopting different Ansatze (14.21) for the effective threshold Seff('r) and hxing these thresholds by minimizing (14.311 : the results are 
presented for the central values of all OPE parameters, (c) The r-dependent effective thresholds as obtained by our algorithm. 
The integer n = 0,1, 2, 3 is the degree of the polynomial in our Ansatz (gJl for Sefr('r): dotted lines (red )—n = 0; solid lines 
(green )—n — 1; dashed lines (blue )—n = 2; dot-dashed lines (black )—n = 3. 



3.25 3.5 3.75 4 4.25 4.5 4.75 5 


fi [GeV\ 



3.25 3.5 3.75 4 4.25 4.5 4.75 5 


yu [GeV] 


(a) (b) 

Fig. 4: Renormalization-scale dependence of the predicted decay constants: (a) /s“'‘'(7i) and /fl*‘‘*(p)i (d) /ss'^'(m) and /^“‘^'(/r). 
For each decay constant, we depict the /r-related uncertainty, i.e., the standard deviation calculated assuming a flat p distribution 
in the range /j, = 3-5 GeV. Dotted lines (red)—vector beauty mesons; solid lines (blue)—pseudoscalar beauty mesons. 


( 99 ). {aGG)) = (192.6 ± S^yst) (^1 - (^1 + 

As is obvious from (14.41) and (I4.7L the OPE uncertainties cancel out, to a great extent, in the ratio, which, consequently, 
can be predicted with a rather high accuracy: 


/b- //b = 0.944 ± O.OlloPE ± O.OlSsyst. 


(4.8) 























Fig. 5: Distributions of the ratios /b* //s and /s* //s^ of beauty-meson decay constants, obtained by generating 1000 bootstrap 
events. For both ratios, their final distributions possess Gaussian-like shapes, with the standard deviations quoted in the plots. 


The main contribution to the OPE error in the ratio arises from the gluon condensate, which enters with different sign 
in the pseudoscalar and the vector correlator (in detail: ±0.01(aGG)i0-005mi,±0.001(qq)). The total uncertainty of the 
ratio is dominated by the systematic uncertainties of the decay constants. Figure [5] shows the distribution of the ratio 
as obtained by a bootstrap analysis. 


B. Decay constant of the meson 


For B *, we choose the same Borel-parameter window as for B* and again require that the deviation of the dual mass 
from the known B* mass does not exceed 10 MeV in the full t window. Our findings for the B*-meson decay constant 
may be cast in the form 


f^'i^\fi = -p,mb,{ss),{aGG)) = (213.6 ±6) ( 1- 


13.2 

213.6'^ 


. 11.8 , 

1 H- S, 

213.6 


{ss) 




MeV, (4.9) 


where /i is defined in (14.111) and 


I (ss) 11/3-0.248 GeV 
0.033 GeV 


(4.10) 


Unfortunately, the sensitivity of fs* to the choice of the scale /i at which the vector correlator is evaluated turns out 
to be rather pronounced. This dependence on the choice of /i may be parametrized by a series in powers of log(/i/7l): 

= 213.6 MeV [1-0.12 log(^/7l)-b 0.11 log^(/r/7l)-b 0.43 log3(^/7l)] , /I = 3.86 GeV. (4.11) 


Averaging over the OPE parameters (using Gaussian distributions of all OPE parameters except for /i, for which a flat 
distribution in the range yL = 3-5 GeV is assumed) yields 

/s. = (213.6 ± 18.2ope ± 6,yst) MeV, (4.12) 

with the following main contributions to the OPE error: 11.5 MeV from the s-quark condensate and 14.1 MeV from Wf,; 
an uncertainty of 3.2 MeV arises from the dependence of Jb*- 
For the pseudoscalar Bg meson, our corresponding estimates read 

(ss), (aGG)) = (225.6 ± 3syst) (l - (^1 + 


MeV. (4.13) 
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As seen in Fig.|4l the sensitivity of to the choice of /x is negligible. The total OPE uncertainty of is rather large: 

/s, = (225.6 ± 18.3ope ± 3syst) MeV. (4.14) 

The decomposition of the OPE error reads: 11.5 MeV are due to the error of s-quark condensate and 14.1 MeV due 
to the error of mb{mb), the uncertainties of the other OPE parameters contribute at the level of 1 MeV. 

Similar to the fs* / fs case, except for the gluon-condensate contribution the OPE uncertainties cancel, to a great 
extent, in the ratio of the decay constants, which may thus be predicted rather accurately: 

= 0.947 ± 0.023ope ± 0.020syst. (4.15) 

The OPE uncertainty in the ratio is dominated by the sensitivity of /b* to the choice of the scale /x. The (obligatory) 
bootstrap analysis gives for the ratio fs* /fsa the nearly Gaussian distribution shown in Fig. [5] 


5. SUMMARY AND CONCLUSIONS 

Exploiting the tools offered by QCD sum rules, we analyzed in great detail the decay constants of the beauty vector 
mesons, paying special attention to the uncertainties arising in our predictions for the decay constants: the OPE error, 
related to the precision with which the QCD parameters are known, and the systematic error, intrinsic to the QCD 
sum-rule approach as a whole, reflecting the limited accuracy of the extraction procedure. Our findings are as follows: 

(i) As was already noted in the case of heavy pseudoscalar mesons Q , also for the vector correlator the perturbative 
expansion in terms of the heavy-quark pole-mass does not exhibit good convergence. Reorganizing the OPE in 
terms of the corresponding running mass allows us to choose a range of scales for which, upon evaluation of the 
correlator, the perturbative hierarchy becomes explicit. Eor scales /x < 2.5-3 GeV, also the running-mass OPE 
does not exhibit any hierarchy of perturbative contributions; at too large scales /x > 5-6 GeV, we observe a strong 
cancellation between the large positive zero-order and the large negative first-order contributions, thus signalling 
that the accuracy of the OPE may deteriorate. There is, however, a sizeable interval of scales, 3 < /x (GeV) < 5, 
where the 0(a^)-truncated OPE provides a good description of the dual correlation function. 

(ii) Requiring the known value of the meson mass to be well reproduced in a relatively broad r window leads, in the 
case of the vector mesons, to some correlation between the scale /x at which the correlator is evaluated and the 
upper boundary of the r window: for /x > 5 GeV, the Borel window for the vector correlator shrinks and thus no 
meaningful extraction of the decay constants of B* and B* from sum rules is possible. The observed correlation 
between the parameters of the Borel window and the value of strongly reduces the (unphysical) /x dependence 
of the extracted beauty-meson decay constants. 

(hi) The r-dependence of the effective threshold and the details of the algorithm for fixing this quantity are crucial 
for obtaining realistic estimates of the systematic uncertainty of the extracted decay constant. For the analysis of 
the ratios of the decay constants of vector to pseudoscalar beauty mesons, where the mass splitting between the 
vector and the pseudoscalar partners amounts to some 45 MeV only, the stringent requirement to reproduce this 
splitting and the individual masses of vector and pseudoscalar beauty mesons with an accuracy not worse than 5 
MeV in the full r range is crucial for obtaining the low systematic uncertainty of the extracted decay constants. 

(iv) The decay constants of pseudoscalar and vector beauty mesons exhibit a strong dependence on the precise value 
oirribijfib). Therefore, the R(s) and B*^-^ decay constants suffer from large OPE uncertainties. The systematic 
uncertainties of the extracted decay constants are of the level of a few MeV and remain under good control. 

(v) The ratios /s* //s and Jb* /can be predicted with very good accuracy because of large cancellations between 
the OPE uncertainties in the ratios and a good control over the systematic uncertainties of the decay constants. 
Our final results read 

/b-//b = 0.944 ±0.021, /bj//b. = 0.947 ±0.030, 

where the error given is the total uncertainty, including the systematic and the OPE uncertainty. The resulting 
distributions are close to normal distributions (Fig.[S]), thus the quoted errors are Gaussian standard deviations. 

(vi) Our results are in excellent agreement with and have a precision comparable to the recent lattice QCD values [l^ 

/b. //b = 0.941 ± 0.026, /bj //b, = 0.953 ± 0.023. 
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